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Abstract 

Time-series  measurements  of  velocity,  temperature  and  conductivity  on  the  northern  California 
shelf  during  two  winter  seasons  permit  an  observational  test,  in  vertically  integrated  form,  of  a 
simple  set  of  subinertial  momentum  and  heat  balances  for  the  bottom  boundary  layer,  which  have 
resulted  from  recent  theoretical  work.  These  are:  (1)  an  along-isobath  momentum  equation  that 
reduces  to  a  classical  Ekman  balance;  (2)  a  cross-isobath  momentum  equation  in  which  the  Ekman 
balance  is  modified  by  a  buoyancy  force  caused  by  distortion  of  the  isopycnal  surfaces  within 
the  boundary  layer;  and  (3)  a  heat  balance  in  which  variability  of  temperature  is  produced  by 
cross-isobath  advection.  The  measurements  confirm  the  importance  of  buoyancy  in  the  cross¬ 
isobath  momentum  equation,  and,  as  has  recently  been  predicted  theoretically,  they  indicate  that 
buoyancy  is  a  dominant  effect  when  the  boimdary  layer  is  thick,  which  typically  occurs  during 
downweUiag-favorable  flows.  An  Ekman  balance  describes  subinertial  fluctuations  in  the  along- 
isobath  momentum  equation  with  only  moderate  success.  In  contrast  to  idealizations  made  in  most 
theoretical  work,  a  buoyancy  force  caused  by  an  along-isobath  temperature  gradient  is  as  important 
as  bottom  stress  in  the  mean  along-isobath  momentum  equation,  and  along-isobath  advection  is 
as  important  as  cross-isobath  advection  in  the  heat  balance. 
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1  Introduction 

Recent  theoretical  work  on  the  bottom  boundary  layer  has  focused  on  a  problem  in  which  a  subin- 
ertial,  along-isobath  flow  begins  from  rest  in  a  stably  stratified  ocean  above  a  gently  sloping,  non- 
conductiug  sea  floor  (see  the  review  by  Garrett  et  al.  [1993]  and  the  recent  articles  by  Ramsden 
[1995a,  1995b]  and  Middleton  and  Ramsden  [1996]).  In  this  problem,  the  sloping  bottom  permits 
the  cross-isobath  T'.kTna.u  transport  to  modify  the  density  field,  and  the  combination  of  the  sloping 
bottom  and  the  distortion  of  isopycnal  surfaces  by  mixing  and  advection  produces  cross-isobath 
buoyancy  forces  that  influence  the  velocity  field.  The  interaction  between  density  and  velocity  leads 
to  evolution  and  structure  that  are  profoundly  different  from  the  more  classical  results  for  a  flat 
bottom.  In  particular,  if  the  fordng  is  steady,  the  models  predict  evolution  toward  a  steady  state 
in  which  the  along-isobath  velocity  is  brought  to  zero  at  the  bottom  by  a  frictionless  thermal-wind 
balance  that  occurs  within  the  boundary  layer,  so  that  the  bottom  stress  and  cross-isobath  Ekman 
transport  vanish  and  the  bottom  offers  no  frictional  resistance  to  the  overlying  flow. 

Although  idealized,  the  recent  theoretical  work  provides  a  concise  set  of  vertically  integrated, 
subinertial  momentum  and  heat  balances  that  can,  at  least  in  principle,  be  tested  observationaUy. 
These  are:  (1)  an  along-isobath  Ekman  balance,  in  which  bottom  stress  is  proportional  to  cross- 
isobath  transport;  (2)  a  cross-isobath  momentum  equation  in  which  the  Ekman  balance  is  modified 
by  a  buoyancy  force;  and  (3)  a  heat  balance  in  which  cross-isobath  advection  produces  temporal 
variability  of  temperature.  These  balances  represent  a  simple  but  oceanographically  relevant  set  of 
boundary  layer  processes,  and  they  are  fundamental  to  the  existing  understanding  of  the  bottom 
boundary  layer  and  its  response  to  the  overlying  flow.  In  particular,  the  concept  of  an  along-isobath 
Ekman  balance  is  a  cornerstone  of  geophysical  fluid  dynamics  (e.g.,  Pedlosky,  1979). 

The  above  balances  have  not,  to  our  knowledge,  been  tested  observationaUy.  In  fact,  few  oceanic 
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studies  have  resolved  the  structure  of  the  bottom  boundary  layer.  Several  researchers  have  measured 
velocity  and  temperature  within  a  few  meters  of  the  bottom,  without  corresponding  measurements 
in  the  upper  part  of  the  boundary  layer  (see,  for  example,  the  reviews  by  Grant  and  Madsen  [1986], 
Dyer  and  Soulsby  [1988]  and  Cacchione  and  Drake  [1990]).  Other  workers  have  reported  coarsely 
spaced  measurements  in  the  upper  part  of  the  boundary  layer,  without  corresponding  measure¬ 
ments  near  the  bottom  (Kundu,  1976;  D’Asaro,  1982a;  D’Asaro,  1982b;  Lentz  and  Trowbridge, 
1991).  Observations  with  profiling  instrumentation  have  resolved  the  entire  boundary  layer,  but 
the  durations  have  been  too  short  for  a  definitive  examination  of  the  sub-inertial  dynamics  (Mer¬ 
cado  and  van  Leer,  1976;  Weatherly  and  Martin,  1978;  Armi  and  D’Asaro,  1980;  Dickey  and  van 
Leer,  1984).  Recent  observations  on  continental  slopes  have  focused  on  the  structure  up  to  hundreds 
of  meters  above  bottom,  and  they  have  not  resolved  the  much  thinner  frictional  boundary  layer  in 
which  Ekman  dynamics  are  relevant  (Thorpe,  1987;  Thorpe  et  al.,  1990;  White,  1994).  Weatherly’s 
(1972)  moored  array  resolved  the  bottom  boundary  layer  in  the  Florida  current,  but  the  duration 
(one  week)  was  too  short  for  a  statistically  significant  examiaation  of  subinertial  variability,  and 
the  analysis  was  limited  primarily  to  estimating  bottom  stress,  determining  the  thickness  of  the 
logarithmic  layer,  and  documenting  Ekman  veering.  Bird  et  al.  (1982)  reported  measurements  that 
spanned  the  bottom  boimdary  layer  on  the  Bermuda  Rise  for  a  period  of  eight  months,  but  the 
vertical  resolution  was  coarse  (four  sensors  between  0.8  and  62  m  above  bottom)  and  the  analysis 
was  limited  primarily  to  a  qualitative  comparison  with  simulations  based  on  a  turbulence  closure 
model. 

We  report  measurements  obtained  as  part  of  the  Sediment  Transport  Events  on  Shelves  and 
Slopes  (STRESS)  program,  which  resolve  the  vertical  structure  of  the  bottom  boundary  layer  on  the 
northern  California  shelf  during  each  of  two  winter  seasons.  The  purpose  of  the  analysis  is  to  test 
the  vertically  integrated  momentum  and  heat  balances  that  have  resulted  from  recent  theoretical 
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work.  In  contrast  to  previous  observations,  the  STRESS  measurements  have  a  duration  sufficiently 
long  to  examine  subinertial  variability  and  they  span  the  boundary  layer  with  vertical  resolution 
sufficiently  fine  to  estimate  the  terms  in  the  vertically  integrated  balances.  In  the  following,  we 
first  present  the  theoretical  background  (Section  2)  and  describe  the  measurements  and  analysis 
(Section  3).  We  then  present  and  discuss  the  results  (Sections  4  and  5)  and  summarize  conclusions 
(Section  6). 

2  Background 

2.1  Mathematical  model 

The  mathematical  model  addresses  an  idealized  configuration  in  which  a  turbulent  boundary  layer 
above  a  non-conducting,  gently  sloping  sea  floor  is  overlain  by  a  stratified  interior  in  which  turbulent 
fluxes  of  momentum  and  heat  are  negligible  (Figure  1).  The  coordinate  system  is  defined  so  that 
the  x-y  plane  coincides  with  the  sloping  sea  floor,  with  x  cross-isobath,  y  along-isobath,  and  z 
perpendicular  to  the  bottom,  with  z  =  0  at  the  bottom.  Time  is  denoted  t.  The  along-isobath 
and  cross-isobath  spatial  scales  are  assrimed  to  be  sufficiently  large  that  advective  accelerations 
are  negligible  in  comparison  with  Coriolis  accelerations.  The  relationship  between  density  p  and 
temperature  T  is  assumed  to  be  approximately  linear;  i.e.,  p  ci  po  -  P{T  -  To),  where  po  is  a  fixed 
reference  density  and  I3  and  To  are  constants.  The  temperature  gradients  dTfdx  and  dT/dy  are 
assumed  to  be  approximately  independent  of  z  within  the  boundary  layer.  This  configuration  is 
essentially  the  same  as  in  the  theoretical  work  reviewed  by  Garrett  et  al.  (1993),  although  it  is 
slightly  more  general  in  that  local  accelerations  are  not  neglected  and  dT/dy  is  not  assumed  to  be 
zero. 

The  dynamics  of  the  boundary  layer  are  described  by  simplified  momentum  and  heat  equations 
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for  a  Boussinesq  fluid  (e.g.,  Pedlosky,  1979).  Tlie  x  and  y  momentum  equations  reduce  approxi¬ 


mately  to 


du  .  dp  ,  dr^ 

Po-qI  ~  Pof'o  =  —ST  +  “al —  ^9P  > 


dt 


dx  dz 


(1) 


and 


dV  _  dp  dXjgy 


The  z  momentum  equation  is  approximately  hydrostatic: 


(2) 


0  = 


(3) 


and  the  heat  equation  reduces  approjdmately  to 


dT  dT  ,  dT  d<j> 
dt^^dx  ^  dy  dz’ 


(4) 


Here  g  is  gravitational  acceleration,  /  is  the  Coriolis  parameter,  a  is  the  small  bottom  slope,  (u,  v,  w) 
is  the  velocity  vector,  p  is  pressure,  Tzx  aJid  r^y  are  components  of  the  stress  tensor,  ^  represents 
turbulent  flux  of  temperature,  and  T,  p  and  po  are  respectively  temperature,  density  and  a  flxed 
reference  density,  a^  before. 

In  (4),  the  term  w  dT/dz  has  been  neglected  without  obvious  justification.  In  the  theoretical 
work  reviewed  by  Garrett  et  al.  (1993),  this  term  is  identically  zero  because,  by  assumption,  the 
velocity  is  everywhere  parallel  to  the  bottom,  so  that  tu  =  0.  Here  we  point  out  that  even  if  w 
is  not  precisely  zero,,  w  dT/dz  is  small  in  comparison  with  udT/dx  within  the  boundary  layer, 
provided  that  the  temperature  in  the  boundary  layer  is  relatively  well  mixed.  To  see  this,  note 
that  w  is  at  most  0{au),  except  in  front-like  regions  with  short  cross-isobath  or  along-isobath 
scales,  which  are  excluded  from  the  analysis.  In  addition,  dT/dx  is  0{adToo/dz),  where  Too  is  the 
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interior  temperature  (e.g.,  Garrett  et  al.,  1993).  Thus  u&T/dx  is  Oiau&Tc^ldz)  while  wdT/dz  is 
0{oLudTldz),  so  that  w  dTIdz  <  udTjdx  provided  that  dTIdz  <  dTcoIdz.  Above  the  boundary 
layer,  where  the  temperature  is  not  well  mixed,  neglect  olw&Tldzis  not  necessarily  justified. 

2.2  Vertically  integrated  momentum  balances 

It  is  convenient  to  define  2:  =  ^  to  be  a  position  just  above  the  turbulent  boundary  layer,  where 
turbulent  fluxes  of  heat  and  momentum  are  negligible.  The  quantity  6  varies  with  time  as  the 
turbulent  boundary  layer  thickens  and  thins. 

To  obtain  a  set  of  vertically  integrated  balances,  begin  by  differentiating  (1)  and  (2)  with  respect 
to  z,  substituting  the  linear  relationship  between  p  and  T,  and  usiug  (3)  to  eliminate  p.  The  result 
is  a  pair  of  equations  describing  the  z  derivatives  of  u  and  u: 


d^u  dv 


-pg 


fdT  dT\  ^ 

Vdx  ^dz)  dz^ 


(5) 


d'^V  ^du  (a.\ 

The  terms  involving  Coriolis  acceleration  and  temperature  in  (5)  and  (6)  represent  the  thermal 
wind  balance,  expressed  iu  a  coordinate  system  in  which  the  x  axis  is  iacliued  at  a  small  angle  with 
respect  to  horizontal. 

Next,  integrate  (5)  and  (6)  from  an  arbitrary  value  of  z  to  z  =  d,  where  dvzxldz  and  dr^yldz 
axe  zero.  The  resialting  expressions  axe 


!d%  rduX  1  „  s  a  j  >  J. 
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and 


where  subscript  6  denotes  evaluation  at  2:  =  ^  and  z'  is  a  dummy  variable  of  integration.  Equations 
(7)  and  (8)  are  a  restatement  of  (1)  and  (2),  with  quantities  involving  p  and  p  expressed  in  terms 
of  us,  vs  and  T.  In  the  quantities  {du/dt)s  and  (dv/dt)s,  the  operations  of  evaluation  at  z  =  5  and 
differentiation  with  respect  to  time  do  not  commute,  because  6  varies  with  time. 

By  integrating  (7)  and  (8)  from  2  =  0  to  z  =  ^,  one  obtains  a  pair  of  vertically  integrated 
momentum  balances: 

and 

where  tsx  and  ny  axe  the  x  and  y  components  of  the  bottom  stress.  The  terms  involving  Coriolis 
acceleration  and  the  stress  terms  in  (9)  and  (10)  represent  the  classical  Ekman  balances.  The  terms 
involving  temperature  represent  buoyancy  forces  that  result  from  non-horizontal  isopycnal  surfaces. 
In  the  special  case  in  which  dT/dy  is  zero  and  local  accelerations  are  negligible,  (9)  and  (10)  are 
consistent  with  the  dynamics  in  the  theoretical  work  reviewed  by  Garrett  et  al.  (1993),  and  they 
axe  identical  to  the  expressions  used  in  the  vertically  integrated  model  developed  by  Trowbridge 
and  Lentz  (1991). 

2.3  Vertically  integrated  heat  balance 

A  simple  form  of  the  vertically  integrated  heat  balance  can  be  obtained  by  integrating  (4)  from 
z  =  0  to  z  =  ^  and  applying  the  condition  =  0  at  the  non-conducting  sea  floor  and  in  the 
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frictionless  interior.  The  result  is 


(11) 


where  we  have  used  the  fact  that  dTjdx  and  dT/dy  are  assumed  to  be  independent  of  z.  The  term 
on  the  left  side  of  (11)  is  the  vertically  integrated  temporal  variability  of  temperature  within  the 
boundary  layer.  The  terms  on  the  right  side  represent  the  effects  of  cross-isobath  and  along-isobath 
advection,  respectively.  Equation  (11)  is  consistent  with  the  theoretical  work  reviewed  by  Garrett 
et  al.  (1993)  in  the  special  case  in  which  dTJdy  is  zero. 

A  more  d3mamically  relevant  form  of  the  heat  balance  involves  the  temperature  anomaly  T—Tg^ 
because  this  quantity  appears  in  the  buoyancy  force  in  (9)  and  because  it  is  related  to  stratification, 
which  influences  vertical  nodxing.  Because  a  heat  balance  involving  T  —  Tg  is  substantially  more 
complicated  than  (11)  and  requires  special  treatment,  we  postpone  consideration  of  this  balance 
until  Section  5. 


3  Methods 

3.1  Measurements 

The  STRESS  observations  (Fredericks  et  al.,  1993;  Trowbridge  and  Nowell,  1994)  occurred  on 
the  northern  California  shelf  (Figure  2)  in  approximately  the  same  location  as  the  Coastal  Ocean 
Dynamics  Experiment  (CODE;  see  Beardsley  and  Lentz,  1987).  The  shelf  in  the  STRESS  region 
consists  of  a  narrow,  steeply  sloping  inner  portion,  between  the  coast  and  the  70-m  isobath,  and 
a  broader,  gently  sloping  outer  portion,  between  the  70-m  isobath  and  the  shelf  break,  at  a  depth 
of  approximately  150  m.  The  primary  site,  denoted  C3,  was  at  a  depth  of  approximately  90  m.  A 
secondary  site  on  the  same  isobath,  denoted  C3',  was  displaced  from  C3  by  approximately  5  km 
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to  the  southeast.  Au  additional  secondary  site,  denoted  C4,  was  at  a  depth  of  approximately  130 
m,  approximately  8  km  olFshore  of  C3.  We  use  a  coordinate  system  in  which  the  y  axis  is  positive 
toward  a  bearing  of  317°  relative  to  true  north,  which  coincides  with  the  isobath  orientation  at 
C3,  and  the  x  axis  is  positive  toward  a  bearing  of  47°.  The  topography  near  C3  has  well  defined 
cross-isobath  and  along-isobath  directions,  and  our  results  are  not  sensitive  to  small  changes  in  the 
orientation  of  the  coordinate  system. 

During  1988-89  (STRESS-l),  A.  J.  Williams,  of  the  Woods  Hole  Oceanographic  Institution, 
deployed  a  bottom  tripod  at  C3,  and  B.  Butman,  of  the  U.  S.  Geological  Survey,  deployed  subsurface 
moorings  at  C3  and  C3'  (Figure  3).  The  tripod  supported  six  BASS  acoustic  current  meters 
and  eight  temperature  sensors  between  heights  of  0.2  and  6  m  above  bottom.  The  subsurface 
mooring  at  C3  supported  five  vector  avera^g  current  meters  (VACMs),  which  record  velocity 
and  temperature,  at  heights  between  6  and  30  m,  and  the  subsurface  mooring  at  C3'  supported 
VACMs  at  6  and  18  m  above  bottom.  Three  of  the  VACMs  at  C3  were  fitted  with  temperature- 
conductivity  sensors  and  optical  transmissometers,  and  both  VACMs  at  C3^  were  fitted  with  optical 
transmissometers.  The  bottom  tripod  weis  within  300  m  of  the  subsurface  mooring  at  C3.  The 
tripod  and  moorings  were  on  site  from  6  December  1988  to  24  January  1989  (STRESS— la)  and 
from  26  January  1989  to  27  February  1989  (STRESS-lb).  All  of  the  instrumentation  functioned 
satisfactorily  during  both  of  these  periods. 

STRESS-1  coincided  with  the  Surface  Mixed  Layer  Experiment  (SMILE;  see  Alessi  et  al.  [1991] 
and  Dever  and  Lentz  [1994]).  As  paxt  of  SMILE,  S.  J.  Lentz  and  R.  C.  Beardsley,  of  the  Woods 
Hole  Oceanographic  Institution,  deployed  a  surface  mooring  at  C3.  This  mooring  supported  vector 
measuring  current  meters  (VMCMs),  which  record  velocity  and  temperature,  between  a  depth  of 
50  m  and  the  water  surface.  We  use  the  SMILE  measurements  briefly  (Section  4.1)  to  illustrate 
the  qualitative  structure  above  the  STRESS  array. 
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During  1990-91  (STRESS-2),  A.  J.  Williams  deployed  bottom  tripods  and  B.  Butman  deployed 
subsurface  moorings  at  both  C3  and  C4  (Figure  3).  At  each  site,  the  tripod  was  within  600  m 
of  the  mooring.  Each  tripod  supported  BASS  current  meters  and  temperature  sensors  between 
heights  of  0.1  and  6  m  above  bottom.  The  mooring  at  C3  supported  an  array  of  VACMs  fitted 
with  temperature-conductivity  sensors  and  transmissometers,  as  well  as  an  array  of  temperature- 
conductivity  sensors.  The  mooring  at  C4  supported  a  VACM,  an  array  of  VACMs  fitted  with 
temperature-conductivity  sensors  and  transmissometers,  a  temperature-conductivity  sensor,  and 
an  array  of  temperature  sensors.  Both  tripods  and  both  moorings  were  on  site  and  successfully 
collecting  data  from  12  January  1991  to  8  March  1991.  Two  of  the  BASS  sensors  at  each  site,  one  of 
the  VACM  velocity  sensors  at  C3,  and  one  of  the  VACM  conductivity  sensors  at  C3  malfunctioned 
(Figure  3). 

All  of  the  STRESS  instruments  obtained  measurements  essentially  continuously  but  recorded 
at  a  number  of  different  rates.  In  the  analysis,  we  use  measurements  that  have  been  averaged  and 
resampled  to  produce  time  series  of  hourly  averaged  values. 

3.2  Estimation  of  the  terms  in  the  vertically  integrated  balances 

The  STRESS  measurements  provide  hourly  averaged  time  series  of  u,  v  and  T  with  relatively 
fine  vertical  resolution  at  C3  and  coarser  resolution  at  C3^  and  C4.  To  estimate  the  terms  in 
the  vertically  integrated  balances,  we  use  the  measurements  at  C3  to  evaluate  quantities  involving 
vertical  structure,  and  we  use  the  measurements  at  C3'  and  C4  to  provide  information  about  cross¬ 
isobath  and  along-isobath  gradients.  To  focus  on  sub-inertial  dynamics,  we  first  estimate  the  terms 
in  the  vertically  integrated  balances  by  using  the  hourly  averaged  data,  and  we  then  apply  filter  pl64 
(e.g.,  Limeburner,  1985),  which  leaves  the  energetic  subinertial  fluctuations  nearly  unattenuated, 
while  removing  the  energetic  tidal  and  near-inertial  fluctuations  (Section  4.1). 
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Some  aspects  of  this  procedure  are  straightforward.  Appropriate  values  of  po,  f  and  a  are  1000 
kg  m“^,  9.1  X  10"®  s“^,  and  0.005,  respectively  (Figure  2).  To  determine  ^  we  use  the  best-fit 
values  obtained  from  the  temperature-conductivity  measurements,  which  are  0.35  kg  m“^  deg“^  C 
for  STRESS-1  and  0.29  kg  m"^  deg"^  C  for  STRESS-2  (Figures  4a  and  4b).  To  put  the  velocity 
and  temperature  measurements  at  C3  on  a  common  grid  suitable  for  evaulating  vertical  integrals, 
we  use  the  heights  of  the  functioning  current  sensors  as  nodes.  For  nodes  at  which  temperature 
measurements  were  not  obtained,  we  assign  a  temperature  record  by  interpolating  between  the 
nearest  adjacent  temperature  measurements.  This  operation  results  in  a  grid  of  eleven  nodes  for 
STRESS-1  and  eight  nodes  for  STRESS-2.  We  use  a  midpoint  rule  to  estimate  integrals  with 
respect  to  z  and  centered  differences  to  estimate  time  derivatives. 

To  determine  S,  we  assume  that  the  turbulent  boundary  layer  coincides  with  the  bottom  mixed 
layer.  We  therefore  set  S  equal  to  the  height  of  the  lowest  sensor  above  the  mixed  layer,  except 
in  the  relatively  rare  occasions  when  the  mixed  layer  extends  to  the  top  of  the  STRESS  array, 
when  we  set  S  equal  to  the  height  of  the  top-most  sensor.  We  estimate  mixed-layer  thickness 
by  using  the  procedure  described  by  Lentz  and  Trowbridge  (1991),  in  which  the  mixed  layer  is 
defined  to  extend  to  the  highest  sensor  at  which  the  local  temperature  is  within  AT  of  the  bottom¬ 
most  temperature,  where  AT  is  an  estimate  of  the  uncertainty  in  the  temperature  measurements. 
Here  we  use  AT  =  0.05°  C,  as  in  Lentz  and  Trowbridge  (1991).  For  the  most  part,  estimates  of 
mixed-layer  height  are  below  the  height  of  the  top-most  sensor  (Section  4.2,  Figures  9  and  10),  so 
that  the  measurements  resolve  the  entire  mixed  layer.  Exceptions  occur  primarily  during  the  first 
half  of  STRESS-2  (Section  4.2,  Figure  10),  when  the  mixed  layer  reached  the  top  sensor  during  a 
number  of  events,  so  that  the  top  of  the  mixed  layer  may  in  fact  have  been  above  the  top  of  the 
measurement  array.  In  Section  4.1,  we  present  a  crude  evaluation,  based  on  estimates  of  gradient 
Richardson  number,  of  the  assumption  that  mixed  and  boundary  layers  coincide. 
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Direct  estimation  of  dT/dx  is  possible  for  STRESS— 2  but  not  for  STRESS— 1.  To  obtain  a  time 
series  of  dT/dx  during  STRESS-2,  we  use  tie  mean  of  tie  estimates  at  nine  ieigits  between  0.19 
and  19  m  based  on  temperature  differences  between  C3  and  C4.  Tie  resulting  time  series  of  dT / dx 
ias  a  mean  of  0.087  ®C  per  km  and  a  standard  deviation  of  0.045  per  km,  so  tiat  fluctuations 
are  typically  smaller  tian  tie  mean,  witi  cianges  in  sign  occurring  less  tian  2%  of  tie  time. 
Tie  time  series  of  dT/dx  at  tie  different  ieigits  are  remarkably  consistent  (temporal  means  and 
standard  deviations  vary  by  less  tian  10%  at  tie  different  ieigits  and  tie  correlation  coeffident  r 
is  greater  tian  0.8  for  all  pairs  of  estimates).  Tiis  consistency  justifies  tie  assumption  (Section  2) 
tiat  dT/dx  is  independent  of  z.  To  estimate  dT/dx  during  STRESS-1,  we  use  a  fixed  value  equal 
to  tie  mean  of  tie  time  series  obtained  during  STRESS-2.  Altiougi  crude,  this  procedure  captures 
tie  dominant  feature  of  tie  bottom  temperatTire  field,  which  is  tie  offshore  cooling  associated  witi 
increasing  depth.  STRESS-1  results  axe  insensitive  to  tie  precise  value  of  dT/dx  tiat  we  use,  and 
STRESS-2  results  are  degraded  but  still  provide  useful  information  (Section  4.2)  if  we  use  a  fixed 
value  of  dT/dx  rather  tian  the  measured  time  series. 

Direct  estimation  of  dT/dy  is  possible  for  STRESS-1  but  not  for  STRESS-2.  To  determine 
dT/dy  during  STRESS-1,  we  use  tie  mean  of  tie  estimates  at  two  ieigits  (6  and  18  m)  based  on 
temperature  differences  between  C3  and  C3'.  Tie  resulting  time  series  of  dT/dy  ias  a  mean  of 
-0.015  ®C  per  km  and  a  standard  deviation  of  0.020  ®C  per  km,  so  tiat  fluctuations  are  comparable 
to  tie  mean.  Tie  time  series  of  dT/dy  at  tie  two  heights  are  only  moderately  weE  correlated 
(correlation  coefiicient  r  =  0.61),  but  they  have  tie  same  dominant  features,  which  axe  a  roughly 
constant  negative  value  witi  brief  periods  of  positive  values. 

To  estimate  bottom  stress,  we  use  tie  log-profile  method  (e.g.,  Sternberg,  1968),  which  is  based 
on  tie  standard  unstratified  wall-layer  model  (e.g.,  Monin  and  Yaglom,  1971),  in  which  tie  flow 
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and  the  flow  coincides  in  direction  with  the  bottom  stress.  Here  u*  is  the  shear  vdocity,  k  is  von 
Karman’s  empirical  constant,  approximately  equal  to  0.40,  and  zq  is  the  bottom  roughness  length. 
To  apply  the  method,  one  fits  velocity  measurements  to  (12),  using  and  (u,/K)ln(zo)  as  free 
parameters,  and  then  determines  the  magnitude  of  the  bottom  stress  from  the  definition  |  |=  po^w*- 

Here  we  fit  vdocity  measurements  obtained  by  the  bottom  three  functioning  BASS  sensors  to  (12), 
and  we  assume  that  the  bottom  stress  coinddes  in  direction  with  the  velodty  averaged  over  the 
bottom  three  sensors.  We  use  three  sensors  because  the  logarithmic  modd  applies  only  at  small 
hdghts  and  because  three  is  the  TuiTiiTmim  number  that  permits  an  assessment  of  the  accuracy 
of  the  stress  estimates  based  on  discrepandes  in  the  modd  fits.  In  addition,  the  lowest  complex 
empirical  orthogonal  function  (EOF;  see,  for  example  Davis  [1976];  here  based  on  the  complex 
vdocity  u  +  iv),  which  accounts  for  more  than  85%  of  the  variance  in  the  hourly  velocity  data  from 
the  entire  current  meter  array,  indicates  that  while  the  speed  profiles  are  logarithmic  throughout  the 
BASS  array,  the  direction  is  vertically  uniform  (particularly  during  STRESS-1)  only  at  the  lowest 
three  sensors,  with  subtle  but  consistent  changes  in  direction  at  the  fourth  sensor  and  above.  Thus 
the  model  (12),  which  implies  a  logarithmic  speed  profile  and  a  vertically  uniform  direction,  appears 
to  apply  only  at  the  bottom  three  sensors.  The  rdationship  between  shear  velocity  and  flow  speed 
that  results  from  the  log-profile  estimates  is  roughly  consistent  with  a  quadratic  drag  law  based  on 
zq  =  0.01  m  (Figures  4c  and  4d),  which  corresponds  to  a  drag  coefficient  c<i  =[  fj  |  I[pq{v?  v^)]  of 
7.5  X  10“®  at  z  =  1  m. 


Log-profile  estimates  of  bottom  stress  are  uncertain  because  of  uncertainties  in  current  mea- 
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surements  and  because  of  sbortcomings  in  the  logarithmic  model  caused,  for  example,  by  distortions 
of  the  near-bottom  velocity  produced  by  acceleration  (Soulsby  and  Dyer,  1981)  and  stratification 
due  to  temperature  or  suspended  sediment  (e.g.,  Monin  and  Yaglom,  1971).  To  quantify  the  effect 
of  errors  in  current  measurements  in  a  simple  way,  assume  for  the  moment  that  the  model  (12) 
is  correct,  that  errors  in  current  measurements  have  a  Gaussian  probability  distribution,  and  that 
the  standard  deviation  of  measurement  errors  is  equal  to  the  standard  deviation  of  the  ensemble 
of  all  residuals  from  fits  of  (12)  to  the  hourly  averaged  velocity  measurements  from  the  bottom 
three  BASS  sensors,  which  is  0.0014  m/s  for  STRESS-1  and  0.0030  m/s  for  STRESS-2.  Given 
these  assumptions,  the  95%  confidence  intervals  for  estimates  of  u,  are  ±0.001  m/s  for  STRESS-1 
and  ±0.002  m/s  for  STRESS-2,  indicating  uncertainties  in  stress  estimates  that  are  much  smaller 
than  the  stresses  themselves  (Figures  4c  and  4d).  Our  opinion  is  that  uncertainties  in  estimates 
are  in  fact  substantially  larger  than  ±0.001  or  ±0.002  m/s,  and  that  they  are  caused  primarily 
by  shortcomings  in  the  logarithmic  model.  However,  calculations  of  the  effect  of  flow  acceleration 
based  on  the  model  of  Soulsby  and  Dyer  (1981),  the  effect  of  stratification  due  to  temperature  based 
on  an  approach  siTini1a.r  to  that  recently  proposed  by  Friedrichs  and  Wright  (1997),  and  the  effect  of 
stratification  due  to  suspended  sediment  based  on  the  model  developed  by  Glenn  and  Grant  (1987) 
indicate  that  these  effects  were,  for  the  most  part,  appreciable  only  when  the  bottom  stress  was 
weak,  so  that  they  are  not  important  for  the  present  purposes.  More  importantly,  with  a  single 
exception  (Section  5.1),  failures  to  close  the  vertically  integrated  momentum  balances  did  not  occur 
during  periods  in  which  any  of  these  effects  appeared  to  be  large. 
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4  Results 

4.1  Qualitative  temporal  variability  and  vertical  structure 

Before  proceeding  to  the  vertically  integrated  balances,  a  brief  examination  of  temporal  variability 
and  vertical  structure  provides  a  useful  qualitative  characterization  of  the  observations.  Here  we 
present  rotary  spectra,  filtered  time  series,  vertical  profiles  of  statistics,  and  estimates  of  gradient 
Richardson  number. 

Rotary  spectra  indicate  that  velocity  fluctuations  occurred  predominantly  in  a  subinertial  band 
(absolute  value  of  frequency  less  than  approximately  0.5  cpd),  a  diurnal  band,  a  near-inertial 
band,  and  a  semidiurnal  band  (Figures  5a  and  5c).  Diurnal  and  semidiurnal  fluctuations  had  pre¬ 
dominantly  clockwise  polarization,  correspondiug  to  negative  frequencies  in  Figure  5,  while  subin¬ 
ertial  currents  were  primarily  along— isobath  and  therefore  had  no  particular  polarization.  During 
STRESS-1,  approximately  55%  of  the  velocity  variance  at  the  top-most  sensor  was  subinertial,  and 
approximately  35%  occurred  ia  the  combination  of  diurnal,  semidiurnal  and  near-inertial  bands.  In 
contrast,  approximately  85%  of  the  variance  at  the  top-most  sensor  during  STRESS-2  was  subin¬ 
ertial.  During  both  STRESS-1  and  STRESS-2,  fluctuations  in  bottom  stress  were  predominantly 
subinertial  (Figures  5c  and  5d). 

Time  series  of  low-pass-filtered  (i.e.,  subinertial)  along-isobath  velocity  indicate  energetic  fluc¬ 
tuations  with  tune  scales  of  several  days,  during  which  the  flow  was  predominantly  poleward  with 
brief  periods  of  equatorward  flow  (Figures  6a  and  6c).  Time  series  of  filtered  temperature  (Figures 
6b  and  6d)  indicate  long  periods  with  strong  stratification  and  other  periods  when  temperature 
throughout  the  measurement  array  was  well  mixed.  The  qualitative  variability  of  temperature  is 
generally  consistent  with  a  simple  model  based  on  the  first  two  terms  in  (11),  in  which  poleward 
flow,  corresponding  to  downwelling  Ekman  transport,  produces  warming,  and  equatorward  flow. 
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corresponding  to  upwelling  Ekman  transport,  produces  rapid  cooling.  However,  there  are  obvious 
exceptions  to  this  simple  view  (e.g.,  day  40  diiring  STRESS-1).  Thick  mixed  layers  generally  corre¬ 
spond  to  periods  of  strong  poleward  flow,  as  reported  by  Lentz  and  Trowbridge  (1991)  for  summer 
CODE  observations. 

Statistics  of  velocity  and  temperature  provide  a  useful  means  of  illustrating  vertical  structure. 
The  mean  cross— isobath  velocity  for  both  STRESS— 1  and  STRESS— 2  is  offshore  near  the  bottom 
and  onshore  above  approximately  20  m  above  bottom,  with  a  uniform  structure  at  heights  above  30 
m  (Figure  7a).  The  vertical  structure  of  the  mean  cross-isobath  velocity  is  remarkably  well  defined 
by  the  measurements,  espedaUy  considering  the  small  magnitude  and  the  fact  that  three  different 
sensors  (BASS,  VACM  and  VMCM)  produced  the  data.  The  mean  along-isobath  flow  is  poleward 
at  all  sensors,  with  a  magnitude  generally  larger  than  that  of  the  mean  cross-isobath  flow  (Figure 
7c).  The  standard  deviation  of  low-pass-filtered  cross-isobath  velocity  (Figure  7b)  is  comparable 
to  the  corresponding  mean,  and  the  standard  deviation  of  low-pass-filtered  along-isobath  velocity 
(Figure  7d)  is  roughly  twice  the  corresponding  mean.  Mean  temperature  profiles  indicate  similar 
stable  stratification  dming  STRESS-1  and  STRESS-2,  although  temperatures  were  larger  during 
STRESS-2  (Figure  7e),  and  standard  deviations  of  low-pass-filtered  temperature  have  little  vertical 
structure  (Figure  7f). 

Calculations  of  gradient  Richardson  number  Ri,  based  on  the  hourly  velocity  and  temperature 
measurements  and  defined  by  Ri  =  Pg{dT/dz)f[{du/dz)‘^+{dvldz)%  provide  a  qualitative  indica¬ 
tion  of  the  likely  importance  of  stable  stratification  in  suppressing  vertical  mixing.  Histograms  of  Ri 
based  on  finite  difference  calculations  for  both  STRESS-1  and  STRESS-2  indicate  values  generally 
well  below  the  critical  value  of  1/4  indicated  by  linear  theory  for  onset  of  instability  (Miles,  1961; 
Howard,  1961)  within  the  bottom  two  meters,  and  values  generally  greater  than  1/4,  with  rare 
excursions  below  1/4,  at  heights  above  roughly  ten  meters  (Figure  8).  Between  heights  of  one  and 
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ten  meters,  the  modal  value  of  Ri  gradually  increases,  and  above  fifteen  meters  the  distribution  of 
Ri  is  approximately  independent  of  height.  These  calculations  support  the  use  of  the  unstratified 
wall  layer  model  (12)  to  estimate  bottom  stress  based  on  measurements  within  roughly  two  meters 
of  the  bottom.  The  calculations  of  Ri  indicate  that  stable  stratification  probably  had  a  significant 
effect  on  vertical  Tui-xiug  at  heights  greater  than  a  few  meters  above  bottom. 

To  evaluate  the  assumption  (Section  3.2)  that  the  mixed  layer  provides  a  reasonable  estimate  of 
the  thickness  6  of  the  turbulent  boundary  layer,  we  estimated  Ri  for  all  adjacent  pairs  of  sensors  by 
differencing  hourly  velodty  and  temperature  measurements,  and  we  defined  an  alternative  boundary 
layer  tliirVTipgs  S^,  as  the  height  of  the  lowest  sensor  at  which  Ri  reaches  unity.  Low-pass-filtered 
time  series  of  6  and  agree  well  in  magnitude  and  axe  weU  correlated  (sample  correlation  coefidcient 
r  =  0.75),  which  indicates  that  the  mixed  layer  does  in  fact  provide  at  least  a  crude  estimate  of 
a  transition  between  a  boundary  layer,  where  shear  instability  and  turbulence  can  occur,  and  a 
non-turbulent  overlying  flow. 

4.2  Vertically  integrated  balances 

To  examine  the  vertically  integrated  balances,  we  regard  the  terms  on  the  right  sides  of  (9),  (10) 
and  (11)  as  forcing,  and  we  regard  the  terms  on  the  left  sides  as  response.  Table  1  summarizes  the 
means  and  standard  deviations  of  the  various  terms  in  these  equations. 

We  first  consider  the  cross-isobath  momentum  equation.  For  both  STRESS-l  and  STRESS-2, 
—pof  Jq  (v—vs)  dz  and  —Ux^  response  and  forcing  in  a  simple  Ekman  balance,  differ  substantially 
iu  magnitude  (Table  1)  and  axe  not  significantly  correlated  at  the  95%  confidence  level  (Figure  9a 
and  10a).  Lack  of  correlation  is  caused  primarily  by  large  fluctuations  in  —pof  JqV  —  vs)  dz,  which 
are  not  matched  by  — Thus  the  measurements  do  not  support  a  conventional  Ekman  balance 
in  the  cross-isobath  momentum  equation.  However,  if  we  include  the  other  terms  in  (9),  the  mean 
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balance  closes  to  within  roughly  35%  (Table  1),  and  the  fluctuations  in  forcing  and  response  are 
comparable  in  magnitude  and  easily  significantly  correlated  at  the  95%  confidence  level  (Figures  9b 
and  10b).  The  dominant  terms  in  the  balance  are  -po  /„  (u  -  vs)  dz  and  Pg{S'^l2)(dT / dx),  and  the 
other  terms  play  smaller  but  non-negligible  roles  (Table  1).  Figures  9b  and  10b  indicate  that  most 
of  the  discrepancies  between  forcing  and  response  occur  at  relatively  high  frequencies  (time  scales 
of  a  few  days  or  less),  and  the  correlations  in  fact  improve  to  r  =  0.87  for  Figure  9b  and  r  =  0.89 
for  Figure  10b,  with  excellent  agreement  in  magnitude,  if  we  apply  a  four-day  boxcar  filter. 

For  both  STRESS-1  and  STRESS-2,  the  dominant  buoyancy  term  in  the  cross-isobath  mo¬ 
mentum  equation,  ^g(S^f2){dT/dx),  is  large  when  the  mixed  layer  is  thick  (Figures  9  and  10). 
In  STRESS-1,  the  buoyancy  force  is  much  larger  than  the  along-isobath  bottom  stress  during 
two  brief  downweUing  events  (i.e.,  events  with  offshore  bottom  stress)  centered  on  days  -12  and 
40  (Figure  9b).  During  STRESS-2,  the  buoyancy  force  is  much  larger  than  the  bottom  stress 
during  two  downweUing  events  of  longer  duration,  centered  on  days  25  and  35,  and  also  during  a 
brief  upweUing  event  (i.e.,  an  event  with  onshore  bottom  stress)  centered  on  day  50.  For  the  most 
part,  these  observations  are  consistent  with  theoretical  work  (e.g.,  Trowbridge  and  Lentz,  1991) 
indicating  that  mixed  layers  are  thick  during  downweUing-favorable  flows,  and  that  a  thick  mixed 
layer  produces  a  large  cross-isobath  buoyancy  force  that  becomes  the  dominant  effect  balancing 
-pQf  joi'o- vs) dz. 

We  next  consider  the  along-isobath  momentum  balance.  In  this  case,  the  terms  in  a  simple 
Ekman  balance,  po/ ^q{v-v.s)  dz  and  -rsy,  dominate  the  fluctuations  (Table  1)  and  are  significantly 
correlated  at  the  95%  confidence  level,  although  they  are  not  weU  correlated  (Figures  11a  and  11c). 
Addition  of  the  terms  involving  temporal  acceleration  and  buoyancy  has  a  small  effect  on  the 
character  of  the  fluctuations  and  on  the  correlation  between  forcing  and  response  (Figure  lib). 
Although  small  in  comparison  with  the  fluctuations,  the  mean  terms  in  (10)  close  remarkable  weU 
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for  STRESS-1,  and  the  mean  along-isobath  buoyancy  force  is  as  important  as  the  mean  along- 
isobath  bottom  stress  (Table  1).  The  presence  of  a  dynamically  significant  buoyancy  force  in  the 
mean  along-isobath  momentum  balance  is  not  consistent  with  the  theoretical  work  reviewed  by 
Garrett  et  al.  (1993),  in  which  along-isobath  variability  was  assumed  to  be  zero. 

Although  statistically  significant  at  the  95%  confidence  level,  the  correlation  between  response 
and  forcing  for  the  heat  balance  (11)  is  small  if  we  consider  only  cross-isobath  advection  (Figures  12a 
and  12c).  The  correlation  between  forcing  and  response  improves  dramatically  if  we  include  along- 
isobath  advection  (Figure  12b),  which  can  be  estimated  only  for  the  STRESS-1  measurements. 
For  STRESS-1,  the  standard  deviations  of  cross-isobath  and  along-isobath  advection  are  nearly 
identical,  and  the  mean  along-isobath  advection  is  substantially  larger  than  the  mean  cross-isobath 
advection  (Table  1).  As  in  the  momentum  balance,  the  importance  of  along-isobath  advection  in 
the  heat  balance  is  not  consistent  with  the  theoretical  work  reviewed  by  Garrett  et  al.  (1993)  in 
which  along-isobath  variability  was  neglected. 

5  Discussion 

5.1  Discrepancies  in  the  momentum  balances 

In  STRESS-1,  discrepandes  in  the  cross-isobath  momentum  balance  occurred  primarily  during 
brief  events  near  day  -10  and  just  before  day  40  (Figure  9)  and  discrepancies  in  the  along-isobath 
momentum  balance  occurred  during  the  same  periods  and  also  near  day  10  (Figure  11).  Am  exam¬ 
ination  of  these  periods  is  instructive. 

During  day  -10,  transmissometer  records  indicate  the  largest  turbidities  that  were  recorded 
during  STRESS-1,  apparently  because  of  wave-induced  sediment  resuspension  during  a  severe 
storm.  Stable  stratification  by  suspended  sediment  might  have  influenced  the  near-bottom  flow 
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during  tHs  period,  resulting  in  over-estimation  of  the  magnitude  of  bottom  stress  by  the  log  profile 
method.  This  idea  is  consistent  with  analyses  by  Gross  et  al.  (1992)  and  Wiberg  et  al.  (1994) 
of  optical  transmission  measurements  obtained  during  the  STRESS  program,  which  suggest  that 
sediment  concentrations  were  large  enough  during  a  few  severe  storms  to  influence  the  near-bottom 
flow.  Also  in  support  of  this  idea,  estimates  of  stress  obtained  from  a  quadratic  drag  law  with  a 
fixed  drag  coefficient  based  on  zq  =  0.01  m  remove,  essentially  completely,  the  discrepancies  during 
this  period  in  both  (9)  and  (10).  Use  of  a  quadratic  drag  law  with  a  fixed  drag  coefficient  does  not 
necessarily  improve  closure  of  (9)  and  (10)  during  other  periods,  probably  because  of  real  variability 
in  the  effective  bottom  roughness,  which  is  captured  by  the  log-profile  method  but  not  by  a  drag 
law  with  a  fixed  drag  coefficient. 

In  the  discrepancy  in  the  along— isobath  momentum  balance  that  occurs  near  day  10  (Figure 
11)^  Po  -  ^s)  dz  is  much  smaller  than  -ny.  This  discrepancy  disappears  essentially  completely 
(without  introducing  a  discrepancy  in  the  cross-isobath  momentum  balance)  if  we  set  6  equal  to 
the  height  of  the  top-most  sensor,  rather  than  the  sensor  just  above  the  mixed  layer.  This  fact 
suggests  that  mixed  layer  thickness  is  a  poor  estimate  of  boundary  layer  thickness  during  this 
particular  event;  i.e.,  that  T^y  just  above  the  mixed  layer  is  not  negligible,  as  required  at  a;  =  5  in 
the  mathematical  model  (Section  2.2).  These  ideas  are  amplified  by  an  examination  of  velocity  and 
temperature  profiles.  The  mean  temperature  profile  during  day  10  (Figure  13b)  indicates  a  mixed 
layer  with  a  height  of  approximately  5  m.  Above  the  mixed  layer,  there  is  a  thick  transition  layer 
between  heights  of  roughly  5  and  20  m  above  bottom,  where  the  stratification  is  much  larger  than 
the  stratification  above  20  m.  The  cross-isobath  velocity  is  nonzero  and  sheared  in  the  transition 
layer  as  well  as  the  mixed  layer  (Figure  13a).  It  is  likely  that  turbulent  momentum  transfer  extended 
above  the  mixed  layer  and  well  into  the  transition  layer  during  this  period,  so  that  mixed— layer 
height  is  a  poor  estimate  of  boundary  layer  height.  We  believe  that  the  vertical  structure  in  Figure 
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13  is  typical  of  strong  upwelling  events,  but  this  hypothesis  is  difficult  to  test  based  on  the  existing 
observations  because  strong  upwelling  events  axe  rare  in  the  STRESS  measurements.  The  vertical 
structure  during  the  more  common  strong  downwelling  events  indicates  a  thick  well-mixed  layer 
overlain  by  a  stratified  interior,  without  a  transition  layer,  so  that  mixed-layer  height  is  probably 
an  accurate  estimate  of  boundary  layer  height  during  downwelling. 

Dever  and  Lentz  (1994)  observed  anomalously  large  horizontal  temperature  gradients  near  day 
40  in  the  SMILE  measurements,  which  had  more  extensive  spatial  coverage  than  the  STRESS 
measurements,  and  they  tentatively  attributed  this  behavior  to  the  presence  of  a  mesoscale  eddy 
on  the  shelf.  It  is  likely  that  during  this  period  the  scales  of  the  density  field  changed  and  that  our 
estimates  of  cross-isobath  and  along-isobath  buoyancy  forces,  which  are  crude  and  involve  implicit 
assumptions  about  spatial  scales,  are  inaccurate.  Failure  to  estimate  buoyancy  forces  accurately 
might  have  caused  the  discrepancies  in  the  momentum  balances  that  occur  near  day  40. 

While  not  definitive,  the  above  discussion  suggests  that  shortcomings  in  our  ability  to  dose 
the  momentum  balances  based  on  the  STRESS-1  measurements  result  from  occasional  problems 
in  estimation  of  bottom  stress,  boundary  layer  thickness,  and  cross-isobath  and  along-isobath 
buoyancy  forces.  During  STRESS-2,  failures  to  close  the  momentum  balances  are  less  event-Hke 
(Figures  10  and  11)  and  not  as  dearly  linked  to  spedfic  explanations,  but  similar  problems  probably 
occurred.  It  is  noteworthy  that  our  results  do  not  suggest  shortcomings  in  the  dynamics  represented 
by  (9)  and  (10),  but  instead  in  our  ability  to  estimate  the  terms  in  these  equations  based  on  the 
existing  measurements. 

5.2  A  balance  for  temperature  anomaly 

The  heat  balance  (9)  is  useful  because  it  permits  us  to  examine  the  relative  importance  of  cross¬ 
isobath  and  along-isobath  advection  in  determining  the  vertically  integrated  variability  of  tern- 
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perature  within  the  boundary  layer.  A  more  dynamicaEy  relevant  balance,  however,  involves  the 
vertically  integrated  anomaly  in  temperature  within  the  boundary  layer,  relative  to  the  temperature 
in  the  overlying  flow.  As  noted  in  Section  2.3,  temperature  anomaly  is  more  relevant  dynamically 
than  temperature  because  the  anomaly  T -Ts  appears  in  the  buoyancy  term  in  the  cross-isobath 
momentum  balance  and  because  it  is  related  to  stratification,  which  influences  vertical  mixing. 


To  obtain  a  vertically  integrated  balance  for  temperature  anomaly  that  is  sufficiently  simple 


for  an  observational  test,  let  h  denote  a  fixed  height  large  enough  that  z  —  h  is  always  above  the 


boundary  layer,  so  that  turbulent  fluxes  of  heat  and  momentum  are  ne^gible  at  z  =  h.  As  in 
Section  2,  we  assume  that  wdT/dz  <  udTIdx  and  that  BT/dx  and  dT/By  are  independent  of  z. 
Equation  (4)  then  yidds 


(13) 


where  subscript  h  denotes  evaluation  sX  z  =  h.  By  subtracting  (13)  from  (4),  integrating  the 
resulting  expression  from  z  =  0  to  z  =  h,  and  applying  the  condition  ^  =  0  at  both  limits,  one 


obtains 

^l\T-n)dz=-§  j\v  -  »*)  d. ,  (14) 

which  is  a  vertically  integrated  balance  for  the  anomaly  T  —  Th.  The  term  on  the  left  side  is  the 
time  derivative  of  the  vertically  integrated  temperature  anomaly,  and  the  terms  on  the  right  side 
represent  the  effects  of  cross-isobath  and  along-isobath  advection.  Note  that  in  (13)  we  can  change 
the  order  of  evaluation  zX  z  —  h  and  differentiation  with  respect  to  t  because  h  is  fixed.  Similarly, 
in  (14),  use  of  fixed  h  permits  changing  the  order  of  integration  with  respect  to  z  and  differentiation 
with  respect  to  t.  Use  of  the  variable  height  z  =  S{t)  would  have  resulted  in  additional  terms  in 
(14),  which  axe  extremely  difficult  to  evaluate  observationally. 

Our  attempts  to  close  (14)  based  on  the  STRESS  measurements,  with  h  fixed  at  the  height 
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of  the  top-most  sensor  or  ajiy  other  sensor,  fail  convincingly.  The  magnitudes  of  our  estimates 
of  the  advective  terms  in  (14)  are  comparable  to  that  of  the  time-derivative  term,  but  the  means 
are  different  and  the  variability  is  essentially  uncorrelated.  The  squared  correlation  coefficient 
between  the  time  derivative  term  and  the  sum  of  the  advective  terms  is  0.01  for  both  STRESS-1 
and  STRESS-2,  which  is  well  below  the  critical  value  required  to  reject  the  null  hypothesis  of 
zero  correlation  at  the  95%  confidence  level.  These  results  do  not  change  if  we  neglect  one  of 
the  advective  terms,  and  they  also  do  not  change  if  we  assume  that  the  subinertial  momentum 
balances  (without  temporal  acceleration)  are  exactly  correct,  so  that  we  can  use  the  right  sides  of 
(9)  and  (10)  (with  6  replaced  by  K)  to  estimate  one  or  both  of  /(f  (u  —  Uh)  dz  and  Sq{v  —  Vh)  dz. 
Failtire  to  dose  (14)  is  probably  a  consequence  of  the  fact  that  the  terms  in  this  balance  involve 
only  differences  m  velodty  and  temperature,  which  are  perhaps  not  accurately  determined  by  the 
measurements,  and  the  fact  that  the  advective  terms  mvolve  dTfdx  and  dTfdy,  which  aie  not  well 
resolved  by  the  measurements.  A  related  issue  is  that  the  assumption  of  vertically  uniform  dTidy 
is  probably  unrealistic,  particularly  above  the  bottom  mixed  layer. 

6  Summary  and  conclusions 

We  have  presented  an  observational  test  of  a  set  of  vertically  integrated,  subinertial  balances  for 
momentum  and  heat  in  the  bottom  boundary  layer,  based  on  time-series  measurements  obtained 
on  the  northern  California  shelf  during  each  of  two  winter  seasons.  As  has  recently  been  predicted 
theoretically,  the  measurements  indicate  clearly  that  the  cross-isobath  momentum  equation  does 
not  reduce  simply  to  a  dassical  Ekman  balance,  but  instead  indudes  a  dynamically  significant 
buoyancy  force  resulting  from  the  distortion  of  isopycnal  surfaces  within  the  boundary  layer.  Also 
as  predicted  theoretically,  the  cross-isobath  buoyancy  force  is  dominant  when  the  boundary  layer 
is  thick,  which  typically  occurs  during  downwelling-favorable  flows.  Subinertial  fluctuations  in 
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the  along-isobath  momentum  balance  axe  crudely  consistent  with  an  Ekman  balance.  The  mean 
along— isobath  momentum  balance  is  significantly  affected  by  a  buoyancy  force  produced  by  an 
along-isobath  temperature  gradient,  and  along-isobath  advection  is  as  important  as  cross-isobath 
advection  in  the  heat  balance. 

Three  conclusions  result  from  our  work.  The  first  is  that  a  buoyancy  force  produced  by  distor¬ 
tion  of  the  isopycnal  surfaces  within  the  boundary  layer  is  a  significant  term  in  the  cross-isobath 
momentum  balance,  and  that  the  buoyancy  force  behaves  in  approximately  the  manner  predicted 
by  the  recent  theoretical  models.  The  second  conclusion  is  that  along-isobath  variability  of  tem¬ 
perature  has  an  important  effect  on  the  bottom  boundary  layer.  The  along-isobath  temperature 
gradient  infiuences  both  the  heat  balance,  where  it  contributes  along-isobath  advection,  and  the 
mean  momentum  balance,  where  it  produces  an  along-isobath  buoyancy  force  that  modifies  the 
cross-isobath  velocity.  Thus  along-isobath  variability  of  temperature  influences  both  the  evolu¬ 
tion  of  the  temperature  field  and  the  structure  of  the  mean  velocity  field.  These  effects  have  been 
neglected  in  the  theoretical  work  on  the  bottom  boundary  layer  that  was  reviewed  by  Garrett  et 
al.  (1993)  and  in  most  numerical  simulations  of  the  entire  shelf  (e.g.,  Allen  et  al.,  1995;  Allen  and 
Newberger,  1996),  although  they  have  been  included  in  the  theoretical  model  developed  recently 
by  Chapman  and  Lentz  (1997).  The  third  conclusion  is  that  even  though  the  STRESS  measure¬ 
ments  provide  longer  duration  and  better  spatial  resolution  than  previous  measurements  in  the 
bottom  boundary  layer,  our  ability  to  close  vertically  integrated  momentum  and  heat  balances 
based  on  these  measurements  is  relatively  poor,  in  the  sense  that  estimates  of  the  forcing  terms 
typically  explain  less  than  50%  of  the  variance  of  the  response  terms.  Although  difficult  to  deter¬ 
mine  definitively,  the  source  of  the  discrepancies  is  likely  inadequacy  in  estimates  of  bottom  stress, 
boundary  layer  thickness  and  horizontal  temperature  and  density  gradients,  rather  than  a  failure 
of  the  simplified  balances  to  represent  the  dynamics  adequately. 
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Table  1:  Statistics  of  the  Terns  in  the  Vertically  Integrated  Balances  (SI  Units) 


j 


STRESS-1 

STRESS-2 

Term 

mean 

std  dev 

mean 

std  dev 

Po  -  {du/di)s] dz 

-0.21  X  10-2 

1.69  X  10-2 

-0.21  X  10-2 

2.25  X  10-2 

-pofSo{v-vs)dz 

+2.79  X  10-2 

3.38  X  10-2 

+3.15  X  10-2 

5.18  X  10-2 

l}g(Sy2){dTfdx) 

+4.09  X  10-2 

3.01  X  10-2 

+4.09  X  10-2 

4.82  X  10-2 

aPgS^iT-Ts)dz 

-1.30  X  10-2 

0.66  X  10-2 

-1.06  X  10-2 

1.73  X  10-2 

'n>x 

+1.00  X  10-2 

1.98  X  10-2 

+1.49  X  10-2 

2.49  X  10-2 

Po  foli^'o/dt)  -  (dv/dt)s]  dz 

-0.41  X  10-2 

1.70  X  10-2 

-0.16  X  10-2 

1.82  X  10-2 

Pof  So  (u- us)  dz 

-1.14  X  10-2 

2.57  X  10-2 

-2.09  X  10-2 

3.21  X  10-2 

l3g{6y2)idTldy) 

-0.74  X  10-2 

0.82  X  10-2 

... 

... 

-0.64  X  10-2 

3.05  X  10-2 

-1.24  X  10-2 

4.90  X  10-2 

f^(dTldt)dz 

+0.90  X  10-® 

3.01  X  10-5 

+0.90  X  10-5 

4.69  X  10-5 

-(dT/dx)j^udz 

+0.49  X  10-5 

2.58  X  10-5 

+0.71  X  10-5 

3.15  X  10-5 

-(dTldy)£vdz 

+1.12  X  10-5 

2.39  X  10-5 

... 

... 
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Figure  Captions 

Figure  1.  Configuration  addressed  in  tlie  tiieoretical  development.  A  non-conducting  sea  floor 
is  inclined  at  a  locally  uniform  angle  a  with  respect  to  horizontal.  A  turbulent  boundary  layer 
adjacent  to  the  sea  floor  is  overlain  by  a  stably  stratified  interior  in  which  turbulent  fluxes  of 
momentum  and  heat  are  ne^^ble. 

Figure  2.  STRESS  study  area  on  the  northern  California  shelf.  C3  is  the  central  site  and  C2,  C3' 
and  C4  are  secondary  sites.  C3  and  C3'  are  at  water  depths  of  approximately  90  m  and  C4  is  at  a 
water  depth  of  approximately  130  m.  The  coordinate  system  is  oriented  so  that  x  is  cross-isobath 

and  y  is  along— isobath  at  C3. 

Figure  3.  STRESS  instrumentation.  During  STRESS-1,  the  heights  of  the  BASS  sensors  were 

0.21,  0.51, 1.10, 1.96, 2.56  and  4.96  m;  the  hdghts  of  the  temperature  sensors  were  0.37,  0.62, 1.10, 

2.01,  2.51,  3.63,  5.00  and  5.84  m;  and  the  heights  of  the  VACMs  were  6, 12, 18,  24  and  30  m  at  C3 

and  6  and  18  m  at  C3'.  During  STRESS-2,  the  heights  of  the  BASS  sensors  at  C3  and  C4  were 

0.39, 0.74, 1.34, 1.94, 2.54  and  4.94  m;  the  hdghts  of  the  temperature  sensors  were  0.19, 1.08, 1.75, 

2.62,  3.76,  4.48,  5.14  and  5.83  m  at  C3  and  0.15,  2.93,  3.77,  4.53,  5.18,  5.96, 13,  23,  27,  41  and  56 

m  at  C4;  the  heights  of  the  temperature-conductivity  sensors  were  16,  22,  25  and  28  m  at  C3  and 

7  m  at  C4;  and  the  heights  of  the  VACMs  were  7,  10,  13,  19  and  31  m  at  C3  and  10,  19  and  71 

m  at  C4.  During  STRESS-2,  the  third  and  sixth  BASS  above  bottom  at  C3,  the  first  and  third 

BASS  above  bottom  at  C4,  the  second  VACM  velocity  sensor  above  bottom  at  C3,  and  the  lowest 

VACM  conductivity  sensor  at  C3  malfunctioned,  as  noted  in  the  text.  We  do  not  use  the  highest 

STRESS-2  VACM  record  at  C4,  which,  at  71  m  above  bottom,  was  far  above  the  bottom  mixed 

V 

layer.  Sensors  above  35  m  are  not  shown  in  the  figure. 

Figure  4.  Temperature-density  relationships  and  drag  relationships  based  on  hourly  data  at  C3. 
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The  STRESS-1  measurements  in  panel  a  were  obtained  from  temperature-conductivity  sensors  at 
heights  of  6,  18  and  30  m  above  bottom.  The  STRESS-2  measurements  in  panel  b  were  obtained 
from  temperature-conductmty  sensors  at  heights  of  10,  13,  16,  19,  22,  25,  28  and  31  m  above 
bottom.  The  solid  lines  in  panels  a  and  b  represent  least-squares  fits  to  the  temperature-density 
relationship,  which  correspond  to  /?  =  0.35  kg  m~^  deg~^  C  for  STRESS— 1  and  /?  =  0.29  kg  m 
deg-^  C  for  STRESS-2.  Panels  c  and  d  show  the  relationship  between  current  speed,  obtained  from 
velocity  averaged  over  the  bottom  three  BASS  sensors,  and  log-profile  estimates  of  shear  velocity. 
The  solid  lines  in  panels  c  and  d  correspond  to  a  quadratic  drag  law  with  a  bottom  roughness 
parameter  of  zq  =  0.01  m. 

Figure  5.  Rotary  spectra  of  velocity  and  bottom  stress.  Positive  and  negative  frequencies  cor¬ 
respond  to  counter-clockwise  and  clockwise  polarization,  respectively.  The  inertial  frequency  is 
approximately  -1.25  cpd.  We  estimated  rotary  spectra  by  detrending  and  demeaning  each  time 
series,  applying  a  Hanning  window,  computing  the  discrete  Fourier  transform  of  the  data  in  the 
complex  form  (x  component)  -I-  i  x  (y  component),  squaring  the  absolute  value  of  the  Fourier  coef¬ 
ficients,  and  averaging  over  17  adjacent  frequencies  to  obtain  spectral  estimates  with  34  degrees  of 
freedom.  We  obtained  spectra  for  STRESS-1  by  concatenating  the  time  series  for  STRESS-la  and 
STRESS-lb.  In  panels  a  and  c,  the  more  energetic  spectra  correspond  to  velocity  measurements 
higher  in  the  water  column. 

Figure  6.  Time  series  of  low-pass-filtered  velocity  and  temperature. 

Figure  7.  Statistics  of  the  low-pass-filtered  velocity  and  temperature  fields  as  functions  of  height 
above  bottom.  Measurements  above  40  m  are  from  the  SMILE  array. 

Figure  8.  Estimates  of  gradient  Richardson  number  Ri  at  four  heights  based  on  hourly  STRESS- 
1  measurements.  The  gradients  are  calculated  based  on  dififerences  between  measurements  at  two 
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different  heights.  The  heights  listed  for  Ri  axe  the  means  of  the  heights  of  each  pair  of  measurements. 
The  pairs  of  heights  are  18  and  12  m  for  the  top  panel,  12  and  6  m  for  the  second  panel,  4.96  and 
2.56  m  for  the  third  panel,  and  2.56  and  0.21  m  for  the  fourth  panel.  Note  that  the  vertical  spacing 
between  the  pair  of  sensors  is  approximately  2.5  m  for  the  bottom  two  panels  and  6  m  for  the  top 
two  panels.  The  solid  lines  correspond  to  Ri  =  1/4. 

Figure  9.  Terms  in  the  cross-isobath  momentum  balance  (9)  estimated  from  STRESS-1  measure¬ 
ments.  For  comparison,  panel  c  shows  the  low-pass-filtered  thickness  of  the  bottom  boundary  layer. 
A  standard  analysis  based  on  Gaussian  statistics  indicates  that  the  critical  correlation  coefficient  r 
between  fordng  and  response  for  significant  correlation  at  the  95%  confidence  level  is  approximately 
0.42  for  the  balances  in  both  panels  a  and  b. 

Figure  10.  Terms  in  the  cross-isobath  momentum  balance  (9)  estimated  from  STRESS-2  measure¬ 
ments.  For  comparison,  panel  c  shows  the  low-pass-filtered  thickness  of  the  bottom  boundary  layer. 
Critical  values  of  r  for  significant  correlation  at  the  95%  confidence  level  are  0.48  and  0.55  for  the 
balances  in  panels  a  and  b,  respectively. 

Figure  11,  Terms  in  the  along-isobath  momentum  balance  (10).  Critical  values  of  r  for  significant 
correlation  at  the  95%  confidence  level  are  approximately  0.35  for  the  balances  in  all  three  panels. 

Figure  12.  Terms  in  the  heat  balance  (11).  Critical  values  of  r  for  significant  correlation  at  the 
95%  confidence  level  are  approximately  0.25  for  the  balances  iu  all  three  panels. 

Figure  13.  Mean  velocity  and  temperature  profiles  for  upwelling  event  near  day  10  during  STRESS-1 
(January  10, 1989;  see  Figure  6). 
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